Exact quantum master equation via the calculus on path integrals 
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An exact quantum master equation formalism is constructed for the efficient evaluation of quan- 
tum non-Markovian dissipation beyond the weak system-bath interaction regime in the presence of 
time- dependent external field. A novel truncation scheme is further proposed and compared with 
other approaches to close the resulting hierarchically coupled equations of motion. The interplay 
between system-bath interaction strength, non-Markovian property, and required level of hierarchy 
is also demonstrated with the aid of simple spin-boson systems. 



The key quantity in quantum statistical dynamics 
is the reduced density operator, p(t) = ti B pr(t), i.e., 
the partial trace of the total density operator over the 
bath space. It is well known that the exact evolu- 
tion of p(t) can be formulated in terms of path integral 
functional assuming the harmonic linear coupling bath 
model. Its numerical implementation is however much 
tedious in comparison with that of the differential quan- 
tum master equation formulations. Tanimura et ulrLiL 
had constructed exact quantum Fokker-Planck equations 
from path integral formulations based on a Gaussian- 
Markovian bath model. More recently, the construction 
of exact quantum master equations has been achieved 
by exploiting various stochastic approaches to dissipative 
dynamicsi 3 i 4 i 5 i 6 i 7 i 8 

In this article, we generalize the Tanimura-Kubo's 
method 1,2 to arbitrary cases and construct an exact dif- 
ferential formulation of p(t) based on rather simple calcu- 
lus on the path integral functional. Let us start with the 
review of the path integral formulation. The total Hamil- 
tonian assumes Hr(t) = H(t) — J2 a QaF a (t). Here, H(t) 
is the deterministic part governing the coherent motion 
of the reduced system, and the system-bath interaction 
is characterized in terms of the coupling between system 
operators {Q a } and stochastic bath operators {F a (t) = 

p(tQ)pl q ; with 
being the thermal bath density operator, 



e JlBt F a e-' hBt }. For an initial p T (t ) 



P 



cq 



we have [h = 1 and {3 = l/(fcsT) hereafter] 



Pit) 



tr F 



U T (t,t Q ;{F a (t)})p(t )e-P hB 
xUl(t,t ;{F a ,(t)})e^p^ 



(1) 



where £/r(i, t ; {F a (t)}) is the Hilbert-space propaga- 
tor with the total Hamiltonian ifr(£)- Noting that 
e f3h vF a {t)e-P h * = F a (t - i/3), Eq. JTJ can be recast as 



p(t) = (U T (t,t ;{F a (t-if3)})p(to)U{(t,t Q ;{F a ,(t)})). 



Here, ( • ) = tr B ( • p^ 1 ). Denote 

p(t)=U(t,t )p(t ), 



(2) 



(3) 



which defines the dissipative Liouville-space propagator 
U{t,to). Let r = {a, a'} be an arbitrary Liouville-space 
representation for the reduced system. The path inte- 
gral formulation for U(T t ,t;T ,t ) = ((T t \U{t,t )\T )) 



U(T t ,t;T ,t ) = / VTe^T[a,a']e~^ a ^ 



(4) 



Here, S[a] and ^[a^'] denote the action and influ- 
ence functionals, respectively. The latter is given as 
[cf. Eq. ©] 

f[a,a'] = (ex P+ { 8 V f dr Q a [a(T)]F a (T - i[3)} 
X a J *o 

xexp_{-zV/ dTQ a ,[a'{T)]F a ,(T)}). (5) 

For the stochastic bath operators {F a (t)} that satisfy the 
Gaussian statistics, the above equation can be recast in 
terms of the standard Feynman- Vernon form»^ 

r\a,cT\ = exp(-V f dT{Q a [a(T)} -Q [a'(r)]} 

x{Q a [a(r)}-Q:[a'(r)}}), (6) 



with 



Q a [a(t)]=Y / f dTC ab (t-T)Q b [a{r)}. (7) 
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Here, C a b{t — r) = (F a (i)Fb(r)) are the bath correla- 
tion functions, which satisfy the symmetry and detailed- 
balance relations C* b (t) = Cf, a (— t) = C a b{t — 

We are now in the position to derive the exact quan- 
tum master equation via a rather simple calculus on the 
above integral functional formulation. For clarity, we 
shall in the following consider only the single dissipative 
mode case, in which the system-bath interaction assumes 
H'(t) = -QF(t). Consider the time derivative of Eq. (gj 
term by term. Firstly, 

|e*l = -,/, 5 H,(o, a ) e ^l s -,H (i ).e*l (8) 

The last identity defines the notion used in this paper. 
We have also 



We have therefore for p n ; n > [cf. Eqs. (jT3|) and (jlljl ] 



d 



-iS[a'] _ ie -iSl°<'].H(t). 



(9) 



Similarly, 
d 



dt 



T = -Q-{(Qt - Qt)F\ + [(Qt - Q t *)f] -Q, (10) 



with 

Q t = f dTC(t~r)Q[a(r)}. (11) 

Together with Eqs. Q and Q), we obtain 

p(t) = -i[H(t),p{t)]-i[Q, Pl (t)]. (12) 
Here, pi(t) = U 1 (t,t )p(t ), or generally, 

p n (t)=U n {t,t )p(t ), (13) 
with Unit, to) being given in terms of path integral as 

,r t 

U n (r t ,t;T ,t ) = VTe lS ^T n [ ai a']e~ lS ^\ (14) 
and 

T n [a,a'} = (-i) n (Q t -Q t *) n nu,a'}. (15) 
Note that p (t) = p{t). We have further [cf. Eq. (fTTTf) ] 

at 

= {Q-J^+^a, a'] - J" n+ i[a, a']-Q} 

-in{C(0)Q-^ w _i[a,o / ]-C*(0)^ B _i[a > a / ]-0} 
+nJ- n [a, a'], (16) 

where 

F n [a,a'] = -iT n ^i[a,a'] \ dr{C(t - t)Q[ol{t)] 

J t 

-C*(t-r)Q[a'(r)}}. (17) 



Pn{t) 



-i[H(t),p n (t)]-i[Q,p n+1 (t)] 

-in [C(0)Qp n ^(t) - C*(0)A,_i(t)Q] 

+np n (t). (18) 



Here. 



p n (t) = U n (t,t )p(t ), 



(19) 



with U n {t, to) being given in terms of path integral as 



U n (r t ,t;T ,t ) = / VTe lS ^^ n [a,a'}e~^ a l (20) 

JTn 



Both p n [Eq. H21) with Eqs. JT1J and (JT3JI] and p„ 
[Eq. (pi with Eqs. (EOJ and O] are of the (2n) th -order 
in the system-bath interaction as their leading contribu- 
tions. In Eq. i|18|) . Pn(^) depends on p n ±i{t) and p n (t). 
The latter however does not belong to the same hierar- 
chy. To proceed, we shall first convert Eqs. l|18|l into a 
set of hierarchically coupled equations of motion (EOM), 
followed by a proper truncation scheme to close the 
resulting EOM. Note that in some special cases, such 
as the driven Brownian oscillatorsi2iiiii£ and the pure- 
dephasing (non-demolishing) dynamics, the exact pi(t) 
can be obtained in terms of po(t) = p(t); thus the non- 
hierarchical exact quantum master equation can be con- 
structed. 

The hierarchicalization of Eqs. I|18fl can be carried out 
with certain forms of bath correlation function. Con- 
sider, for example, a Gaussian-Markovian bath with 
C(t) — rje~' ft , where 7 is real and r\ may be complex. 
In this case, p n {t) = — JPn(t) and Eqs. i|18[l become com- 
pletely hierarchical, recovering the Fokker-Planck equa- 
tion obtained previously by Tanimura and co-workersi^ 

Consider another model in which the bath spectral 
density assumes J(u>) oc ui/(lu 2 + 7 2 ) 2 , and the result- 
ing bath correlation is found to be of C(t) = fote -7 * + 
j/ie~ 7t , if the Matsubara terms can be neglect ed>i2ii^ 
Here, 7 is real and vq and v\ are complex. The hierar- 
chically coupled EOM equivalent to Eq. Q18[l can read- 
ily be constructed as follows. For the present model 
of C(t), we have Q t = V0Q0 + viQi [cf. Eq. (|TT|) ]. 
where Q = J^dr (t - r) e - 7 (*- T )Q[a(r)] and Q 1 = 
f tQ dTe~^ {t - T ^Q[a{T)]. Consequently, T n of Eq. O can 
be expressed in terms of J^f = Q%(Q' ) k ' ' QliQ'xY ' T ', 

with k + j + k' + j' = n. Introduce now p\j simi- 
larly as Eqs. i(T3|) and ((Hfl but with T n there being re- 
placed by . Note that dQo/dt = Q 1 — jQ and 
dQi/dt = Q t - 7Q1. The hierarchical EOM for 
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with Pqq = p, can then be obtained as 

ft? = -i[H(t),p%f}-(k + j + k>+j'hp>°f 

k'j' 



-LO.^iPfej+i -viPkj \ 



/ fc'-ij'+i 



(21) 



Obviously, pjjiy' is of the [2(fc + j + fc' + j')] th -order in 
the system-bath interaction as its leading contribution. 
In general, the bath correlation function C(t) can be ex- 
pressed in terms of a series of exponential expansion, and 
thus, by following the similar procedure as shown above, 
one can construct a set of hierarchically coupled EOM 
based on Eq. (|T%jl . 

To complete the formulation, we propose in the fol- 
lowing a novel truncation scheme and compare it with 
other two existing approaches to close the hierarchical 
EOM equivalent to Eq. I|18|). As shown earlier p n {t) in 
Eq. fl%|l is of the same order as p„(t), we need only to 
approximately anchor pat in terms of p n <N, such that 
the system-bath interactions are accounted rigorously for 
up to (2N) th order, but partially for higher-order con- 
tributions. To do that, let us first recast Eq. I|15(l as 
T n = —i(Qt — Q t *)T n -i, and then approximate the in- 
volving Q t with its bath-free-evolution counterpart; i.e., 



Qtrun(t) = / dTC(t-T)G(t,T)Q. 



(22) 



't 



Here, G(t, r) denotes the dissipation-free propagator that 
satisfies dQ{t,r)/dt = —i£(t)Q(t,T), with the reduced 
system Liouvillian C(t)- = [H(t),-]. The resulting 
truncation scheme proposed here amounts therefore to 
[cf. Eqs. ® and JE|] 



PN{i) 



Qtmn{t)pN-l(t) - PN-l(t)Ql un (t) . (23) 



The above truncation leads to a time-local evolution for 
the anchoring pjv, and the resulting EOM is in a partial 
ordering prescription (POP). 

In comparison, the corresponding chronological or- 
dering prescription (COP) of (2iV) th -order truncation 
can be realized via setting pN+i = 0, as it effectively 
leads to a time-ordered, memory description of pN{t) 
in terms of |O n <jv( T < t). It is easy to verify that the 
conventional second-order memory-kernel and time-local 
formulationsiiiSiiSiii are just the special cases of the 
present COP-truncation (p2 = 0) and POP-truncation 
[Eq. (gSl) with N = 1] , respectively. 

Consider now the truncation scheme proposed by Tan- 
imura and co-workers in which the fast modulation 
or Markovian anzatz is assumed applicable at the an- 
chor level. In this case, Qtrun(i) of Eq. I|22|) reduces to 



Q™r _ qq^ where c = J™(ItC{t). This scheme is 
essentially identical to that proposed by Tanimura and 
co-workers, but is obtained here for general C(t) with- 
out invoking the Gaussian-Markovian bath model. Un- 
like the iV-level POP [Eq. (E3 with Eq. JUJ] or COP 
(pn+i = 0) truncation, the TV-level fast- modulation 
truncation scheme is not a rigorous (27V) th , but rather 
(2iV — 2) th -order formulation. We shall come back to 
this point later in terms of the practical efficiencies of the 
aforementioned three truncation schemes; see discussions 
following Fig.^ Obviously, all these schemes become ex- 
act when N — > oo. 

To investigate the efficiencies of involving truncation 
schemes for the study of non-Markovian dissipation with 
arbitrary system-bath interaction, we consider here a 
simple spin-boson system: H — ^Qa x , with Q = a z 



and C(t) 



-7* 



In the following, we adopt the 



integrated coupling strength T = J^°dtC(t) = A 2 /7 
and the dimcnsionlcss modulation parameter k = 7/ A 
to characterize the nature of system-bath coupling, and 

s N = imax {ti} | J^ +T dt\A N (t) - /Ia^oo^)) 2 j; with 

T = 2n/il, to calibrate the finite TV-truncation induced 
error in evaluation of a given test quantity A(t), which 
will be chosen to be P22(t)—pn(t) for demonstration. The 
initial condition is set to be Pjk(0) = fajSjk', j, k = 1, 2. 

Depicted in Fig.QJare the resulting error indicators sjv 
for the aforementioned three truncation schemes, demon- 
strated as functions of anchoring index N at various val- 
ues of system-bath coupling strength and modulation pa- 
rameter: (a) T = fi, k = 1; (b) T = O, k = 0.1; (c) 
r = 10O, k = 1; and (d) r = 10 O, n — 0.1. This fig- 
ure clearly demonstrates the following features: (i) The 
anchor index depends not only on the system-bath 
coupling strength T, but more importantly on the mod- 
ulation parameter /?; (ii) In the slow modulation (k « 1) 
limit, the low order truncation may not be sufficient even 
in the weak coupling (r <C f2) regime; (iii) As inferred 
from their constructions, these three truncation schemes 
are in principle of the same quality in the fast modulation 
limit; see Fig.^c). Considering further the computation 
efforts involved, while the anchor p^ in the POP or fast- 
modulation truncation is expressed directly by Eq. ill'-il) 
or its Markovian counterpart, respectively, the shall 
however be propagated via coupled EOM. Thus, the POP 
truncation scheme proposed in the work is overall the 
best, including the case studied in Fig.njb). 

Shown in Fig.[2arc the approximate results of the test 
quantity pTi{t) — p\\ (t), obtained via the POP (solid), the 
COP (dash), and the fast-modulation (dot) truncation 
schemes at the specified value of N in each panel such 
that the error tolerance of 10 -3 5 is met by the minimum 
sat among these three truncation schemes. Included in 
each panel of this figure is also the exact result (thick- 
solid) that can be obtained via any truncation scheme 
with a sufficiently large N. Note that Fig.|2[b) is de- 
picted for 26 < t < 27 (in the unit of 2ir/£l) where the 
difference between various curves is relatively large. The 
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overall superiority of POP-truncation scheme is again 
highlighted. We have also showed (iv) The behavior of 
p(t) in short time, long time, and Markovian regimes may 
be accounted for properly by the low order p n ; (v) The 
relatively high order p n are required for the intermediate 
time, and the slow modulation regimes; (vi) Note that 
the critical damp occurs at T = 17 in the fast modula- 
tion (re ^> 1) limit [cf. Fig.^a)], it occurs however at a 
much larger T in the slow modulation (re <C 1) limit [cf. 
FigJJJb)]. 

To summarize, we have derived an exact quantum 
master equation formalism via a simple calculus on the 
Feynman- Vernon influence functional path integrals. It 
is valid for arbitrary external field driven reduced dynam- 
ics under non-Markovian system-bath interaction beyond 
the weak coupling regime. The present derivation consti- 
tutes an alternative approach to the exact quantum mas- 
ter equation, Eq. (jl8(l , which is identical to that obtained 
recently by Shao via rather advanced stochastic differ- 
ential equation algebrai^i^ It is also noticed that Shi 
and Geva had also recently constructed a formally exact 
quantum dissipation theory^ It however depends prac- 
tically on the path-integral evaluation of the Nakajima- 
Zwanzig dissipation kernel, despite the fact that is ex- 
pressed in terms of force-force correlation function to be 
evaluated with correlated system-bath ensemble^ 



We have also proposed a novel truncation scheme; i.e., 
the POP-scheme of Eq. 1123H . which is shown to be over- 
all superior to the two existing truncation approaches. 
This superiority is of much more implication than that 
demonstrated in the numerical examples of this work, as 
in general C(t) contains multi-exponential terms, espe- 
cially at low temperature regime, and the relevant num- 
ber of p n ; with n < N [cf. Eq. i|21|) ] in the hierarchical 
EOM increases quasi-exponentially with the truncation 
anchor N. The numerical demonstrations are also of rich 
physical implications on the interplay between system- 
bath coupling strength, non-Markovian property, and the 
required order of truncation, which had been summa- 
rized in details in Comments (i)-(vi) following Fig.Q]and 
Fig-El Especially, low-order quantum dissipation theo- 
ries should be used with care if non-Markovian dynamics 
is important. 
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re = 1; (b) T = 0, re = 0.1; (c) Y = 10 fi, re = 1; and (d) 
L = 10fi, re = 0.1. 



FIG. 2: The evolutions of pi?(t) — pn(i) implied in Fig.0 
with the specified anchor N in each panel being chosen such 
that the tolerance of 10~ 3,5 is satisfied by the minimum sn 
among the three truncation schemes. The exact evolutions 
can be obtained with any schemes truncated at sufficiently 
high level. The curves for the exact and the POP almost 
overlap with each other in each of the four panels. 
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